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Abstract. We present a model for the hydrogenation 
states of Polycyclic Aromatic Hydrocarbons (PAHs) in 
the diffuse interstellar medium. First, we study the abun- 
dance of hydrogenation and charge states of PAHs due 
to photo-ionization, photo-dissociation in the interstellar 
UV field, electron recombination and chemical reactions 
between PAH cations and H or H 2 . For PAH cations, we 
find that the dehydrogenation effects are dominant. The 
hydrogenation state of PAHs depends strongly on the H 
density, the size of the molecule and UV field. In diffuse 
clouds with low H density and normal UV radiation, PAHs 
containing less than 40 C are completely or strongly dehy- 
drogenated whereas at high H density, they are normally 
hydrogenated. The partially dehydrogenated species dom- 
inate in intermediate density clouds. PAHs above 40 C are 
quite stable and are fully hydrogenated, which would fa- 
vor their spectroscopic search in near IR surveys of Diffuse 
Interstellar Bands (DIBs). 



^ ' Key words: ISM: molecules - clouds - dust, extintion 

■ ■ ' 1. Introduction 

In recent years, it has been suggested that a population of 
large carbonaceous molecules, such as Polycyclic Aromatic 
Hydrocarbons (PAHs) and fullerene compounds could 
play a significant role in the physics and chemistry of inter- 
stellar and circumstellar environments. These molecules 
have been proposed as the sources of the broad infrared 
emission features (Duley & Williams 1981, Leger & Puget 

1984, Allamandola, Tielens & Barker 1985) and of the dif- 
fuse interstellar absorption bands (Leger & d'Hendecourt 

1985, van der Zwet & Allamandola 1985, Crawford, Tie- 
lens & Allamandola 1985). PAHs would represent the most 
abundant molecules in the interstellar medium after H 2 
and CO (Leger & d'Hendecourt 1985). They could exist 
under various charge and hydrogenation states depend- 
ing on local conditions (density, UV flux...). Environment 
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studies of some DIBs suggest carriers with ionization po- 
tential of 10-13 eV such as PAH cations (Cami et al. 1997, 
Sonncntrucker et al. 1997, 1999). 

The aim of this work is to propose a description of the 
abundance of different hydrogenation states for some small 
and medium PAHs < 50C under the physical conditions 
present in the interstellar medium, especially in diffuse 
environments where UV photons strongly influence the 
dissociation of PAHs. 

We first compare the effect of physical processes in- 
volving PAH cations and find the balance dominated by 
photo-dissociation and hydrogenation processes. Then we 
calculate the distribution of hydrogenation states for re- 
presentative PAH cations. 

2. Physical processes of PAHs in diffuse clouds 

We discuss here the processes involving PAHs in diffuse in- 
terstellar clouds. There are four main processes in compe- 
tition: photo-dissociation following the absorption of UV 
photons, photo-ionization, electronic recombination and 
reactions between PAH cations and H or H 2 . 

2.1. Photo-dissociation 

Photo-dissociation takes place when the photon absorbed 
by the molecule has a sufficient energy to dissociate the 
molecule. We only take into account the dissociation of a 
H or H 2 of PAH cations. 

After an UV photon is absorbed, the PAH cation is 
either dissociated, or stabilized by the emission of IR pho- 
tons. The UV photon energy is distributed on all vibra- 
tional degrees of freedom of the molecule. When the ener- 
gy accumulated in a vibration mode becomes comparable 
to the C— H dissociation energy (~ 4 eV for PAH cations), 
the molecule dissociates. 

Photo-dissociation studies are often analysed using 
statistical models like the RRKM (Rice-Ramsperger- 
Kassel-Marcus) theory, where one supposes the excited 
molecule is quickly relaxed into the electronic fondamen- 
tal state by internal conversion. Here, we use the photo- 
dissociation rates calculated by Allain et al. (1996) (cf. 
Table 1) 
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Table 1. Table of photo-dissociation rates of PAH cations 
and neutral PAHs derived from Allain ct al. 1996 and 
addition rates of PAH cations from Scott et al. and Le 
Page et al. 



PAH 



cation 



neutral 





H loss 
s" 1 


R~2 loss 
s" 1 


H loss 
s" 1 


H2 loss 
s" 1 


k a dd,H 

cm 3 s _1 


Benzene 


3.67(-9) 


2.07(-9) 


2.89(-9) 


1.50(-9) 


2.2(-10) 


Anthracene 


5.67(-9) 


2.13(-9) 


3.40(-9) 


1.19(-9) 


1.6(-10) 


Pyrene 


3.67(-9) 


1.47(-9) 


2.05(-9) 


7.13(-10) 


1.4(-10) 


Coronene 


4.33(-9) 


1.37(-9) 


1.49(-9) 


4.55(-10) 


l.l(-lO) 


Ovalene 


7.33(-10) 


2.0(-10) 


1.78(-10) 


4.85(-ll) 


9.6(-ll) 


PAH 50C 


9.0(-16) 


1.4(-16) 


2.31(-16) 


3.61(-17) 


8.8(-ll) 



2.2. Reactions between PAH cations and H or Hi 

According to Le Page (private communication), for small 
PAH cations, the Hydrogen addition rate might be propor- 
tional to that of benzene multiplied by the ratio between 
the number of reactive C and the total number of C. The 
reactive C atoms have one bond with a H. The addition 
rates for some small PAH cations are given in Table 1. 
The reaction rate between an ion and a neutral can be es- 
timated by the Langevin rate (~ 10~ 9 cm 3 s _1 ) which is 
larger than the kinetic reaction rate between two neutrals 
(~ lCT 11 cm 3 s- 1 ). 

2.3. Photo-ionization 

Photo-ionization depends on the cross section of UV pho- 
ton absorption and on the local UV intensity. The photo- 
ionization rate is given by: 

47rF ion(a) ( J E)cT UV abs( Q ) {E)F(E)dE. 

The ionisation yield llon=fion/cuv is given by Ver- 
straete et al. (1990). We take the ionisation potential 
IP(Z,N C ) from Dartois and d'Hendecourt (1997). The UV 
intensity F (E) is given by Gondhalekar et al. (1980). 



2.4- Electronic recombination 

We calculate the recombination rate using the cross sec- 
tion of recombination for plane molecules from Spitzer 
(1978) and Verstraete et al. (1990). We find: 



fc roc (cm 3 s- 1 ) = 1.66 x 10" 



+6.3 x 10-« x x ft. 



Even though the Spitzer's approximation may be one 
order of magnitude lower than the recombination rate of 
the naphtalene cation (mesured to be 3 x 10~ 7 cm 3 s _1 
by Abouelaziz et al. 1993), this theory may be justified 
for relatively large PAHs (above 30-50 C), in absence of 
experimental studies. 



2.5. Respective roles of PAH processes 

Interstellar diffuse cloud environments are characterized 
by low H densities (0.1 — 100 cm~ 3 ) and an important UV 
flux (Gondhalekar, 1983). We suppose the temperature 
constant and ~ 100K. The electronic density is supposed 
to be equal to the C density and all C atoms are supposed 
to be ionized in diffuse clouds with n c /nu = nc/^H — 1-4 
10- 4 . 

We now use these main processes rates to estimate 
the abundance of different species of PAHs under different 
charge and hydrogenation states. We suppose a stationary 
state. The abundances are normalized to the total abun- 
dance: 

EilCmHi] = 1. 

The density of each species is given by the equilibrium 
equation: 



d[C m H+] 



dt 



loss + gains. 



d[C m H+] 



dt 



-(krecfle + fcdiss.H + &diss,H 2 + fcadd.HnH 



+ fcadd,H2 n H 2 )[C m H+] + foon[C m H;] + fcdiss.H [C m j ] + 

fcdiss,H 2 [CmH^j] + fc a dd,HWH[C m H^L j] + & a dd,H 2 «H 2 [C m H^L 2 ] , 

where i represents the hydrogenation state of PAH, 
n e , riH and njj 2 respectively the electron, H and H2 den- 
sity, k rcc electronic recombination rate between PAH + and 
electron, fcdi ss ,H and fcdi ss ,H 2 photo-dissociation rates after 
UV photon absorption with ejecting a H or H2, fc a dd,H 
and A: a dd,H 2 H and H 2 addition rates of PAH cations, fc; on 
ionisation rate of neutral PAH. 

We resolve these linear equations using a IDL (Interac- 
tive Data Language) programme. We find that the PAHs 
tend to be dehydrogenated in physical conditions such as 
in diffuse cloud environments where the Hydrogen density 
is about 0.1 — 100 cm~ 3 , the UV intensity is assumed to 
be constant (Gondhalekar intensity) and the temperature 
is about 100K. In what follows, we will then consider the 
two dominant processes concerning hydrogenation states 
of PAH cations: photo-dissociation and hydrogenation. 

3. Hydrogenation states of PAH cations 

The hydrogenation state of PAHs depends on two fac- 
tors: first, the hydrogenation rate, which is function of 
Hydrogen and PAH cation density and second, the photo- 
dissociation rate which strongly depends on the size of the 
molecules and the UV intensity. The equilibrium state is 
a balance between the photo-dissociation rate due to UV 
photon absorption and the hydrogenation rate. We will 
present a simple model which takes only these two pro- 
cesses into account. 



[C m H+] 



[C m H+_ J ] 
- [C r 



3+] 



- [C m H 
[C 
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C 2 ,H 2 
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C 24 H n 
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Fig. 1. Abundance of hydrogenation states of the 
coronene cation as a function of xjj- Three phases are 
present with increasing the H density. At xh < 0.8, 
strongly dehydrogenated species dominate. At xh around 
1 partially dehydrogenated species dominate. At xh > 1-5, 
the coronene cation is completely hydrogenated. 



cation C 32 H+, the PAH 40C cation and the PAH 50C 
cation. We note the decrease of nn with PAH sizes, except 
for coronene which has a high dissociation rate because of 
its particularly symmetric configuration, which limits ra- 
diative relaxation emission of infra-red vibrational modes. 

3.2. Case 2: variable weighting o/fc a( jd vs H cover 



We assume: 



[C„,H,][<iv] _ fc, 
[C m H,_,][H] fcd 



n-(i-l) 



In this case, each hydrogenation state of PAH cations 
is weighted by a factor 2^ in order to take into account 
the number of sites available for adding a H atom. 

The abundance for each hydrogenation state is given 

by: 

Jg»M = ( XH y x x .... x^lxl. 

The number of molecules normalized by the total 
abundance is: 



j=i n— (j— 1) ' 



= 1. 



EIUIc-h*] = [c m H ] (i + E: =1 (2H) l ir;.z; 

We have calculated the balance between hydrogenation 
and photo-dissociation. The results are given in Fig. 2. 



C m H„ + UV ^fdd C m H„_ j + H 

The goal of the model is to estimate the distribution 
of hydrogenation states for a given PAH cation as a func- 
tion of environment parameters in the interstellar medium 
clouds. 

3.1. Case 1: with constant weighting of fc a dd vs H cover 
We assume for 1 < i < n the balance: 

[C„,H,][$, y ] _ fc add 

where fc a dd and fcdiss represent the hydrogenation rate 
and the photo-dissociation rate, and [H] are the local 
UV flux and the H density. 

In this case, we assume that the probability of captur- 
ing or loosing a H is the same for various hydrogenation 
states of PAHs. Then in the second case, we will assign a 
weight to each hydrogenation state. 

Let us put xa = 1^ x -E. One has: 



The number of molecules normalized by the total 
abundance is: 

Er= [c m H i ] = EIUMic m H ] = i. 

The abundance for each hydrogenation state is: 

[C ro Hi] = (xH) l [C m H ] = (mY ^fuy 
xh = 1 corresponds to the maximum of intermediate 
species for n H = 67, 50, 74, 14, 1.7 x 10" 2 , 2 x 10~ 5 cm" 3 
respectively for the anthracene cation Ci4H^~ , the pyrene 
cation CigH^, the coronene cation C24H^ 2 , the ovalene 



4. Discussion 

4-.1. Dependence on H density and on molecule size 

The hydrogenation state strongly depends on the H den- 
sity. Two processes are competing: that of the UV radia- 
tion which leads to the dehydrogenation of PAH cations 
and the reaction between PAH cations and H which favors 
hydrogenated states. 

As we can see in the figures 2, three phases exist for 
each hydrogenation state of PAHs except for PAHs > 40C. 
The UV intensity is assumed to be constant while the 
H density varies. At low density, PAH cations are com- 
pletely or strongly dehydrogenated but at high density, 
they are normally hydrogenated. At intermediate density, 
PAH cations are partially dehydrogenated. 

For PAH cations > 50C, the completely hydrogenated 
state is favored (cf. Fig.^d). This can be explained by the 
fact that large molecules have more degrees of freedom 
available to store the excess energy brought by the UV 
photon, which decreases the photodissociation rate. 

4-2. Search for PAHs and link with DIBs 

Small PAHs < 30C will be difficult to search because their 
abundance is distributed over a number of hydrogenation 
states. Their superimposed absorptions can contribute to 
the envelope of quasi continuum extinction in addition to 
dust due to a large number of states and isomers. PAHs 
> 40C are stable and fully hydrogenated. Their signa- 
ture can be searched in the fundamental band transition 
expected in the far red and near infra-red around 1 /im 
(Foing & Ehrenfreund 1994). Rotational contours for such 
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Fig. 2. Different hydrogenation states of a) the pyrene cation Ci6H^ , b) the coronene cation C24H ] f 2 , c) the ovalene 
cation C32H^ 4 and d) the PAH 50C cation as a function of H density nn in the presence of UV (Gondhalekar spectrum 
multiplied by a constant factor $„). At low densities, strongly dehydrogenated species dominate whereas at nn > 
20, 30, 5 cm~ 3 respectively for the pyrene cation, the coronene cation, the ovalene cation, partially dehydrogenated 
species are important. The normal hydrogenated PAH 50C cation is stable over a large range of densities. 



PAH > 40C would lead to FWHM broadening of 2-4 
cm -1 depending on the exact size and geometry of the 
molecule (Ehrenfreund & Foing 1996). These transitions 
can be searched in near IR surveys of DIBs. Laboratory 
spectroscopy and dissociation and hydrogenation rates are 
also needed for guiding the search for PAHs in space. 
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